Parallel finite volume simulation of the spherical shell dynamo with pseudo-vacuum magnetic boundary conditions

In this paper, we study the parallel simulation of magnetohydrodynamic (MHD) dynamo in a rapidly rotating spherical shell with pseudo-vacuum magnetic boundary conditions. A second-order finite volume scheme based on collocated quasi-uniform cubed-sphere grid is applied to spatial discretization MHD equations. To ensure solenoidal condition field, adopt widely-used approach whereby pseudo-pressure introduced into induction equation. The temporal integration split by approximate factorization approach, resulting two linear algebraic systems both solved preconditioned Krylov subspace iterative method. multi-level restricted additive Schwarz preconditioner domain decomposition and multigrid method then designed improve efficiency scalability. Accurate numerical solutions benchmark cases are obtained our code, comparable existing local results. Several large-scale tests performed Sunway TaihuLight supercomputer show good strong weak scalabilities noticeable improvement from up 10368 processor cores.